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Abstract 

We measure the absolute branching fractions of D s semileptonic decays where the hadron in the 
final state is one of (j), rj, rf, Kg, K*°, and /o, using 2.8 x 10 5 e + e~ — > D S D* decays collected in 
the CLEO-c detector at a center-of-mass energy close to 4170 MeV. We obtain B{Df — > </>e + f e ) = 
(2.29 ± 0.37 ± 0.11)%, B(Df -> r/e+z/ e ) = (2.48 ± 0.29 ± 0.13)%, B(D+ -» r/e + z/ e ) = (0.91 ± 0.33 ± 
0.05)%, where the first uncertainties are statistical and the second are systematic. We also obtain 
B(D+ -> K°e + u e ) = (0.37 ± 0.10 ± 0.02)%, and i?(Z)+ -> K*°e + u e ) = (0.18 ± 0.07 ± 0.01)%, 
which are the first measurements of Cabibbo suppressed exclusive D s semileptonic decays, and, 
B{D+ -» f e + u e ) x B(/ -» 7T+7T-) = (0.13 ± 0.04 ± 0.01)%. This is the first absolute product 
branching fraction determination for a semileptonic decay including a scalar meson in the final 
state. 



* Deceased 
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The study of D s semileptonic decays provides interesting information on several aspects 
of heavy quark decays. First of all, the total semileptonic width provides discrimination be- 
tween different theoretical evaluations of hadronic matrix elements affecting charm semilep- 
tonic decays. The Operator Product Expansion (OPE) predicts that all the charmed mesons 
have the same semileptonic width, modulo non-factorizable corrections [l[ . The ISGW2 form 
factor model [2[ predicts a difference between the D and D s inclusive rates, as the spectator 
quark masses m u and m s differ on the scale of the daughter quark mass m s in the Cabibbo 
favored semileptonic transition. The D + and D° semileptonic widths are equal within the 
3% accuracy of the measurements, and the compositions of their inclusive spectra are domi- 
nated by the lowest lying resonances [3|. This result is explained by the observation that the 
s quark in the final state is usually produced with a small enough momentum to be bound 
to the spectator anti-quark in an I = sq meson [4]. D s semileptonic decays share these 
kinematic features, and thus an absolute measurement of D s semileptonic decays sheds some 
light also on inclusive processes. Specific decays contribute valuable information on light 
meson properties. For example, the fraction of semileptonic decays going into rj and rf is 
sensitive to the pseudoscalar mixing angle, and may indeed shed some light on rj — r/-glueball 
mixing 0. In addition, decays including nir and KK in the final state can elucidate the 
nature of exotic light scalar mesons jl, 0]. 

D s exclusive semileptonic decays have been studied by ARGUS, CLEO, BaBar, and fixed 
target experiments. No absolute measurements of branching fractions exist. The branching 
fraction — > 0£ + z/£, which is the most widely studied jl, 0, [Io|, EH 12], is generally 



normalized with respect to the decay D s — > (pit. However, the Dalitz plot for this mode 
shows the presence of a significant broad scalar resonance whose contribution to the observed 
yields changes depending upon the selection criteria 13j. For this reason, this mode is not 
suitable for normalization. Recently, the BaBar collaboration 14j used the normalization 
mode D s — > KKw with a mass cut of ± 10 MeV around the nominal as suggested in 
Ref. |l3], to obtain B(D+ -> <pe + v e ) = (2.61 ± 0.03 ± 0.08 ± 0.15)%. CLEO measured Q 
the ratio [T(r]£ + u e ) + r(r]'£ + u e )}/[T((f)£ + u e )} = 1.67 ± 0.17 ± 0.17. Finally, BES reported 
the inclusive branching fraction B(Df — > e + anything) = (7-7-^-2'l)%- The uncertainties 
are too large to allow a meaningful comparison between inclusive and exclusive channels. 

We use a data sample of 310 pb" 1 , collected at a center-of-mass (CM) energy close to 4170 
MeV, with the CLEO-c detector [l3, The momenta and directions of charged particles 
are reconstructed in the tracking system, which also provides charged particle identification 
information based on specific ionization (dE/dx). A Ring Imaging Cherenkov Detector 
(RICH) completes the charged particle identification system [19(, and is critical near 1 
GeV, where the specific ionization bands of the K and 7r overlap. The photon energy and 
direction are measured in the Csl electromagnetic calorimeter, whose energy measurement 
E, combined with the momentum p information from the tracking system, provides the key 
electron identification variable E/p. The Csl calorimeter measures the electron and photon 
energies with an r.m.s. resolution of 2.2% at E = 1 GeV and 5% at E — 100 MeV. 

At Ecm — 4170 MeV, the cross section for e + e~ annihilation into D* + D~ + D+D*~ 
is approximately 0.9 nb, while other charm production totals ~ 7 pb, and the light quark 
continuum cross section is ~ 12 nb [20J. We look for semileptonic decays of the in events 
in which the D~ is fully reconstructed in a hadronic mode (tagged events). Each event also 
must include at least one isolated photon, as either the or D~ originates from a D*. 
Here and throughout the paper, charge conjugate decays are implied. Charged tracks are 
used to form the D~ if their fitted helical trajectory approaches the event origin within a 
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distance of 5 mm in the azimuthal projection and 5 cm in the polar projection (9), where 
the azimuthal projection is in the bend view of the solenoidal magnet. In addition, each 
track must possess at least 50% of the hits expected, must be within the fiducial volume 
of the drift chamber, and must have a momentum of at least 40 MeV. Pions and kaons 
are identified using dE/dx and RICH if their momenta are above 700 MeV, otherwise only 
dE/dx identification is used. We form rj and 7r° candidates from pairs of photons that 
deposit energy in the calorimeter in a manner consistent with an electromagnetic shower 
and are not matched to tracks. We require that the two photons to have less than 3cr pull 
mass which is defined as the standard deviation from the expected tt° or r\ mass. In the best 
calorimeter region (| cos#| < 0.71) we use photons with energies greater than 30 MeV, while 
in the endcap region (0.93 > | cos#| > 0.85) we require an energy greater than 50 MeV. 

The tag modes used are listed in Table I. Some tag specific selection criteria are applied. 
For the modes D~ — > K + K~n~ and Dj — > K + K~tt~7t°, the 7r is required to have a 
momentum greater than 100 MeV to suppress the background from D* decays. Similarly, for 
the mode D~ — > 7r + 7r~7r~, two pions of opposite charge must have momentum greater than 
100 MeV. The Charged track pairs used to reconstruct K s (via K s — > 7r + 7r~) are required 
to have an invariant mass within 3a of the K s mass. In addition, for the D~ — > K S K~ 
and Dj ^K*~K*° tags, we require candidate K s to originate from a vertex displaced from 
the interaction point, and the K s momentum vector, obtained from a kinematic fit of the 
charged n momenta, must point back to the beam spot. For resonance decays we select 
intervals in invariant mass centered on the resonance masses [llj and within ±150 MeV for 
p~ — > n~n°, ± 100 MeV for K* — > Kit, ± 10 MeV for rf — > or rf — > fry. In addition, 

for the latter rf decay mode, we apply a helicity angle cut which is the angle measured in 
the rest frame of the decaying parent particle between the direction of the decay daughter 
and the direction of the grandparent particle as | cos^l < 0.8. Tags are required to have 
momentum consistent with coming from D S D* decay. 

TABLE I: Tagging modes and number of signal and background (Bkg) events, determined from 
two-Gaussian fits to the invariant mass distributions. The signal window is ±2.5cr of the D~ mass 
for all modes, except rjp~ (±2er). The last two columns show the corresponding estimates of the 
number signal and background tags accompanied by a reconstructed 7 from D* — > ^D s transition, 
which are determined from the MM* 2 spectrum (see text). 



Mode 


Invariant Mass 


MM* 2 






Signal 


Bkg 


Signal 


Bkg 


K+K-Ti- 


13952 ± 232 


11280 


8245 ± 245 


13970 


K S K- 


2943 ± 128 


561 


1749 ± 146 


1555 


r/ir 


1806 ± 120 


4747 


1241 ± 123 


3936 


r/(?77r + 7r _ )7r~ 


1231 ± 55 


415 


907 ± 109 


1036 


K-K+-K-1T 


5300 ± 401 


34419 


2913 ± 289 


24985 


7r + 7r~7r~ 


4331 ± 716 


25824 


2439 ± 558 


16619 


K*-K*° 


1565 ± 114 


1442 


841 ± 87 


2440 


rjp- 


4002 ± 254 


22044 


2168 ± 268 


18450 


ir'rj'(frf) 


2515 ± 342 


18593 


1817 ± 212 


12061 


Sum 


37645 ± 978 119325 22320 ± 792 95052 
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We further select tags using the recoiling mass M r , 



M rec = y (Eqm - ^ D7 ) 2 " (Pcm - P Ds -) 2 , (1) 

where E CM (pcm) is the CM energy (momentum), E Da (poj is the tag energy (momentum). 
D* daughter tags peak broadly in M rec due to the presence of the photon in the tag side, 
while directly produced D~ tags have a narrow peak. We accept events for which M rec is 
within -55 MeV < M rec — Mp* < 55 MeV. This cut is broad enough to encompass both 
the narrow peak associated with D s tags and almost all the broad component associated 
with tags that originate from D*. Then we fit the invariant mass distribution Md s of these 
events, using a two-Gaussian shape for the signal plus a polynomial background shape. The 
signal component allows us to define an effective a = /i0"i + (1 — /i)o"2 where o\ and a<i 
are the standard deviations of the two Gaussian components and f\ is the fractional area of 
the first Gaussian. We require that the candidate invariant mass to be within 2.5cr (2a for 
the rjp mode) of the nominal D s mass [21] . Random D s backgrounds are estimated through 
sideband samples. We then combine the tag with a well reconstructed 7 and calculate the 
missing mass squared MM* 2 , the square of the invariant mass of the system recoiling against 
the 7-tag pair. 

MM* 2 = (Eqm - E D - - E,) 2 - (p CM - p D7 - p 7 ) 2 , (2) 

where E 1 (p 7 ) is the energy (momentum) of the additional 7. In order to improve the 
MM* 2 resolution, we use a kinematic fit that constrains the D s decay products to Md s and 
conserves overall momentum and energy. 

Fig. [T] shows the MM* 2 distribution for the nine tags considered. In order to estimate 
the number of tags used for further analysis, we use a two-dimensional binned maximum 
likelihood fit of the measured MM* 2 and Md s distributions. We consider three components in 
the fit: a signal, comprising true tags accompanied by the photon from the D*, a background 
composed by true tags combined with a random photon, and a second background comprising 
false tags. We infer the D s combinatoric background from two 5a (4a for the rjp mode) wide 
intervals on both sides of the Md s signal peak. The MM* 2 signal fit is improved by using a 
probability distribution function (PDF) derived from fully reconstructed D S D* events. We 
fit signal and sideband intervals in Md s simultaneously. The sideband intervals constrain 
the shape of the random D s background. We then extract the tag yield from the fitted 
signal function integrated within ± 2.5a around the MM* 2 most probable value. For the 
nine modes considered, Table [I] shows the number of signal and background tags, as well as 
the signal and background tags reconstructed in conjunction with an isolated photon . 

We next describe reconstruction of the semileptonic decays. For any given tag-photon 
combination, we seek a candidate e + and a set of hadrons. Positrons are identified on 
the basis of a likelihood ratio constructed from three inputs: the ratio between the energy 
deposited in the calorimeter and the momentum measured in the tracking system, the specific 
ionization dE/dx measured in the drift chamber, and RICH information [22|. Our selection 
efficiency averages 0.95 in the momentum region 0.3-1.0 GeV, and 0.71 in the region 0.2-0.3 
GeV. The average fractions of charged ir and K incorrectly identified as positrons averaged 
over the relevant momentum range are approximately 0.1%. We study events containing 0, 
1], r)', Kg, K*°, and /o in the final state. Track and 7 selection criteria, as well as resonance 
cuts are the same as used in the tag reconstruction except that we select candidates with 



invariant masses within ± 10 MeV of the known 6 mass 21 for 6 — > K + K , and within 



± 75 MeV of the known K*° mass for K*° — > K + tt . Among the 77 candidates not used 
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FIG. 1: The MM distribution from events with a photon in addition to the D~ tag. In each plot, 
the D~ mode studied is indicated, and the solid curve represents a fit to a Crystal Ball function 
(signal), the dotted curve represents the total background, and the dashed curve represents the 
background composed of random D s tags, constrained by the sideband sample. Both terms are 
well described by 5 th order Chebychev polynomial background functions. 

in forming a tag, we choose the one with the smallest pull mass to form r]e + v e and r]'e + v e 
candidates. We use 7/ — > 7]7t + tt~ only. For the channel K*° — > K + tt~, the K and e must 
have the same charge. Finally, we form candidates by combining two pions of opposite 
charge and require that their invariant mass be within 100 MeV of the known fo mass. In 
each case we require that the event have no unused tracks, and that the tag and semileptonic 
candidate have opposite charge. 

For each 7 candidate, we perform two kinematic fits, one assuming that the 7 combines 
with the tag to form a D*~ , the other assuming that the semileptonic decay comes from 
a D* + parent. We require the D S D* pair to conserve energy and momentum in the CM 
frame, and the mass of the candidate D s to be consistent with the known mass. When we 
assume the tag to be the daughter of a D*~, we constrain the energy of the photon plus 
tag candidate to be consistent with the expected D*~ energy, otherwise we constrain the 
energy of the tag candidate to be consistent with the D~ energy in the CM system. Finally 
we choose the photon and hypothesis with the smallest x 2 and calculate the missing mass 
squared MM 2 defined as 

MM 2 = (E* CM - E* D - -E;- El - EU 2 (3) 

-(-PdT ~P~<-Pe- Phad) > 

where E* (p*) is the energy (momentum) of the positron candidate and E^ ad (Ph ad ) is the 
energy (momentum) of the hadron candidate in the CM system. For signal events, MM 2 
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is the v e invariant mass squared and thus it peaks at zero. Fig. [2] shows the measured 
MM 2 for each final state summed over all tag modes used. We require signal events to 
have a |MM 2 | < 0.05 GeV 2 . We estimate the background coming from random D~ tags 
by studying the sideband samples, and the remaining background by studying a sample of 
simulated DD events that is 20 times bigger than our data set. Fig. [2] shows the signal and 
background distributions for the six modes studied. The MM 2 shapes are well modeled 
by the signal Monte Carlo simulation. We show also the background estimates from the 
sideband data sample, as well as the background predictions from a Monte Carlo simulation 
of charm meson decays at this center-of-mass energy. Note that the background is small 
in all the modes considered. We have also investigated backgrounds produced by random 
photons associated with a true semileptonic event, and we found these to be even smaller, 
thus we do not subtract them. 
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FIG. 2: The MM 2 distribution for tagged semileptonic events in the exclusive modes: Df — > <fte + v e , 
D+ -> r)e + v e , D+ -» r]'e + u e , D+ K*°e + u e , D+ -» K°e+v e: and D+ f e + v e . The solid entries 
correspond to the events from signal regions, the hatched entries are the events from sideband and 
the dashed entries represent the scaled background events from generic MC. 



We evaluate exclusive branching fractions for semileptonic decays including the hadron i 
through the relationship 



Bi 



t (n a ^ Q! bkg) 



e SL (^a n a)Bi 

where the index a runs over the tag modes, the index % represents a specific hadronic final 



had 



(4) 
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state, Bf ad identifies the branching fraction for that state, and eg L represents the average 
efficiency for finding the exclusive semileptonic decay in the tag sample used. We evaluate the 
semileptonic efficiency by considering two Monte Carlo (MC) samples. The first contains 
a tag event accompanied by a semileptonic decay (double-tag sample), while the latter 
contains a tag accompanied by a generic D s decay (single-tag sample). The ratio of the 
single and double tag efficiencies is the desired esL for each tag. Table [IT] shows the signal, 
and background yields, esL and the branching fractions determined for the six semileptonic 
channels. We derive esL from each tag mode independently, and then we compute their 
average weighted by tag abundance. The efficiency obtained with this method accounts for 
the different tag efficiency in semileptonic and generic D s decays. We treat the exclusive 
channel <pe + v e slightly differently: as the </> reconstruction efficiency is strongly momentum 
dependent, we perform the analysis in five 200 MeV wide momentum bins. We attribute 
all the signal events in the K + K~e + v e in the final state to form the 4>ev e channel. Recent 
BaBar studies [3] have estimated the 5-wave fraction to be (0.221q;os) % of the decay rate, 
smaller than our statistical uncertainty in our (f)ev e branching fraction. We also look for tttt, 
KK, and Kit outside the mass windows corresponding to the resonant states studied and 
we found no evidence for additional channels. 



TABLE II: The signal and background yields, the semileptonic efficiency eg L and the derived 
branching fractions for the six semileptonic channels studied. The Df — > foe + u e branching fraction 
quoted represents the product branching fraction B(Df — > foe + v e ) x B(fo — > ir + ir~), which is the 
dominant decay mode in Ref. 2l|. 



Signal Mode 



n 



n 



Df 
Df 
Df 
Df 
Df 
Df 



SL 



(%) 



B{%) 



<pe+v e 45.50 0.06 17.79 ± 0.33 2.29 ± 0.37 
r)e+v e 82.49 0.32 37.65 ± 0.27 2.48 ± 0.29 
ie + v e 7.50 0.06 21.04 ± 0.22 0.91 ± 0.33 
K°e+u e 13.99 0.29 33.14 ± 0.26 0.37 ± 0.10 
K*°e + v e 7.50 0.18 27.52 ± 0.23 0.18 ± 0.07 
foe + v e 13.99 0.88 46.79 ± 0.31 0.13 ± 0.04 



We consider several sources of systematic uncertainty. The dominant component is asso- 
ciated with the number of tags, which is affected by the lack of our knowledge on the random 
7 in the background PDFs. We estimate it by repeating the fit with a variety of shapes, 
namely polynomials of different order, or special shapes derived from MC simulation, and 
obtain an uncertainty of 3.6%. Systematic uncertainties associated with hadron selection 
such as tracking (0.3% per charged particle), K and n identification (0.6% and 0.3% respec- 
tively), and i] selection criteria (2%) have been studied extensively [23j. Similarly, the K s 
selection criteria are derived from Ref. (24J, and have (0.8%) uncertainty. The systematic 
uncertainty on the electron identification efficiency (1%) is assessed by comparing radiative 
Bhabha samples, Bhabha events embedded in hadronic events, and MC samples. The re- 
quirements that there are no extra tracks in the event and that the net charge is zero have 
been evaluated with a data sample comprised of two hadronic tags. The comparison between 
results obtained with this sample and corresponding MC samples give an overall systematic 
uncertainty of 0.6% from these two requirements. Finally, we consider the dependence of 
the efficiency for semileptonic decays on the form factors. The CLEO MC uses the form 
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factors predicted by the ISGW2 model [3]. We have generated also samples based on simple 
pole form factors and compared the efficiencies derived with the two methods to estimate 
this effect. The related systematic uncertainty ranges from 0.1% to 2.4%. 

We check the normalization of our branching fractions by measuring the well known 
branching fraction B(D° — > K~e + v e ) using a D~ D* + sample from the same data set. We 
reconstruct DD* — > D~ D* + — > D~tt + D° decays, where the D° decays into K~e + u e , and the 
D~ decays into these six hadronic exclusive final states: D~ — > K + tt~tt~ , D~ — > K + tt~7t~7t° , 
D - -> #<j7r~, IT -> KgTT Tt^, D~ -> ^7r-7r + 7r~, D - -> ^+7^7^. The selection criteria 
and analysis procedure are the same as used in reconstructing the -D s semileptonic decays. 
We get in total 14759 ± 203 of tagged events and 350 ± 18 signal events, and using Eq. (j4j), 



70. 



we derive a branching fraction B(D — > /<" e + z/ e ) = (3.45 ± 0.21 

This result is in agreement with the two most recent absolute measurements: B(D° — > 
K~e+u e ) = (3.61±0.05±0.05)% from CLEO-c @, based on 281 pb" 1 data at the ^(3770), 
and B{D° -> R-e + u e ) = (3.45 ± 0.07 ± 0.20)% from Belle @ . 



All the measurements reported here are first absolute measurements of exclusive semilep- 
tonic D s decays, moreover this is the first report of Cabibbo suppressed final states and scalar 
meson above the threshold to decay in the observed final state. For the six D s semileptonic 
decays considered, Table II III shows the derived branching fractions including the systematic 
errors. 

TABLE III: The derived branching fractions including the systematic errors for the six semileptonic 
channels studied. The Df — > foe + v e branching fraction quoted represents the product branching 



fraction B(D+ — > /oe + i^ e ) x B(fo ~^ 7r+7r ) 5 which is the dominant decay mode in Ref. [2l| 



Signal Mode B{%) 

D+ -» (f)e + v e 2.29 ±0.37 ±0.11 
D+^i]e + v e 2.48 ±0.29 ±0.13 
Df -» r/'e+fe 0.91 ± 0.33 ± 0.05 
D+^K°e+u e 0.37 ±0.10 ±0.02 
Df -» K*°e+u e 0.18 ± 0.07 ± 0.01 
Df -» foe + v e 0.13 ±0.04 ±0.01 



These results allow us to draw several interesting conclusions. The sum of the branching 
fractions measured imply B(Df — > Xe + u e ) = (6.47 ± 0.60)%, about 16 % below the value 
7.82±0.13)%, inferred from measured D + and Dz inclusive semileptonic branching fraction 



27 ] and the charmed meson lifetimes j21j . We have searched for additional hadronic final 
states formed with two charged tracks, as well as from two charged tracks and a n°, and found 
no evidence for semileptonic decays including other hadronic final states. No other significant 
branching fraction is expected. This result is consistent with the predictions of the ISGW2 
model 0, supporting the conjecture that SU(3) is broken in charm semileptonic decays. On 
the other hand, the difference in widths may arise from non factorizable contributions at 
the level of ~ 10% Q. The ratio B{Df -> r/e+z/ e ) /£(£>+ -> r/e+z/ e ) = 0.36 ± 0.14, is in 



agreement with the previous CLEO result |15J. The ISGW2 model involves a 77/77' mixing 
angle close to —10°, which is the minimum value obtained from mass formulae [2l| if a 
quadratic approximation is used. According to Ref. [EJ , the measured ratio is consistent with 
a pseudoscalar mixing angle of about —17°, provided that a glueball component probability 
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of the order of 10% is present in the r( . Finally, we have the first measurement of a D s 
semileptonic decay including a scalar meson above the threshold for decay to the observed 
final state, which opens up the exciting possibility of elucidating the nature of exotic light 
mesons jof. 
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